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Abstract 

Background: Thousands of long noncoding RNAs (IncRNAs) have been reported in mammalian genomes. These RNAs 
represent an important subset of pervasive genes involved in a broad range of biological functions. Aberrant expression of 
IncRNAs is associated with many types of cancers. Here, in order to explore the potential IncRNAs involved in hepatocellular 
carcinoma (HCC) oncogenesis, we performed IncRNA gene expression profile analysis in 3 pairs of human HCC and adjacent 
non-tumor (NT) tissues by microarray. 

Methodology: Differentially expressed IncRNAs and mRNAs were detected by human IncRNA microarray containing 33,045 
IncRNAs and 30,215 coding transcripts. Bioinformatic analyses (gene ontology, pathway and network analysis) were applied 
for further study of these differentially expressed mRNAs. By qRT-PCR analysis in nineteen pairs of HCC and adjacent normal 
tissues, we found that eight IncRNAs were aberrantly expressed in HCC compared with adjacent NT tissues, which is 
consistent with microarray data. 

Conclusions: VJe identified 214 IncRNAs and 338 mRNAs abnormally expressed in all three HCC tissues (Fold Change >2.0, 
P<0.05 and FDR <0.05) with the genome-wide IncRNAs and mRNAs expression profile analysis. The IncRNA-mRNA co- 
expression network was constructed, which may be used for predicting target genes of IncRNAs. Furthermore, we 
demonstrated for the first time that BC017743, ENST00000395084, NR_026591, NR_015378 and NR_024284 were up- 
regulated, whereas NR_027151, AK056988 and uc003yqb.l were down-regulated in nineteen pairs of HCC samples 
compared with adjacent NT samples. Expression of seven IncRNAs was significantly correlated to their nearby coding genes. 
In conclusion, our results indicated that the IncRNA expression profile in HCC was significantly changed, and we identified a 
series of new hepatocarcinoma associated IncRNAs. These results provide important insights about the IncRNAs in HCC 
pathogenesis. 
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Introduction 

The incidence of hepatocelluar carcinoma is currently the fifth 
highest in solid tumor and is the third leading cause of cancer 
death worldwide, accounting for approximately one million deaths 
annually [1,2]. Geographic areas with highest frequency are 
located in Africa and Eastern Asia. Most cases of HCC are 
attributed to either chronic viral hepatitis infection (hepatitis B or 
C) or cirrhosis [3]. Unfortunately, the survival rate of HCC 
patients remains poor despite recent advances in medical 
treatment and surgical techniques [4]. More seriously, the 
mortality of HCC is still growing with an increasing trend of 
new occurrence [4]. 

The overwhelming majority of human transcriptome was 
confirmed to be noncoding genes [5]. Over the past decade, 
those abundant transcripts have been declared to have important 



regulatory potential in biological processes. Increasing evidence 
have pointed to the critical regulatory role of noncoding RNAs 
(ncRNAs) in normal cellular physiological processes as well as the 
contribution of aberrant ncRNA expression to cancer biology [6] . 
According to their length, ncRNAs can be divided into two major 
categories, long noncoding RNAs which are tentatively defined as 
a class of RNA molecules longer than 200 nucleotides and short 
noncoding RNAs (18 to 200 nucleotides) respectively [7]. It is well 
documented that microRNAs are aberrantly expressed in many 
types of cancer. MicroRNA-21, microRNA-122 and microRNA- 
657 are involved in the development of HCC [8-11]. Recently, 
IncRNAs have been shown to exert critical roles in a series of 
biological processes, including genetic imprinting, immune 
response, tumorigenesis, cellular development and metabolism 
through comprehensive mechanisms [5,12]. LncRNAs also played 
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critical roles in a variety of human diseases, including cancer, 
Huntington's disease and Alzheimer's disease [13-15]. Deregula- 
tion of IncRNAs has been proposed to be associated with 
hepatocarcinogenesis, such as MVIH, HI 9, HEIH, HULC, 
TUC338 and MEG3 [16-22]. For example, IncRNA high 
expression in HCC (IncRNA-HEIH) facilitates tumor growth 
through Enhancer of Zeste Homolog 2 in humans [18]; LncRNA- 
HULC which is a highly specific up-regulated IncRNA in HCC, 
envisaged as a novel biomarker because it could be detected in 
blood of HCC patients [19]; The expression of LncRNA- 
MALAT-1 up-regulated predicts tumor recurrence of HCC after 
liver transplantation [20]. Despite these exciting development, 
many more IncRNAs playing crucial roles in HCC remain to be 
clarilied. 

By analyzing the global expression profile of IncRNAs and 
mRNAs, we identified 214 IncRNAs and 338 mRNAs that were 
significantiy dilferentially expressed in HCC tissues and adjacent 
NT tissues. LncRNA classification and subgroup analysis, genomic 
location analysis, construction of the IncRNA-mRNA co-expres- 
sion network and qRT-PCR analysis were performed for further 
analyze these differentially expressed IncRNAs (Figure 1). We 
found eight IncRNAs were dysregulated in nineteen pairs of HCC 
samples compared with adjacent NT samples and expression of 
seven IncRNAs was significantly correlated to their nearby coding 
genes. Simultaneously, the IncRNA-mRNA co-expression network 
was constructed, which may be used for predicting target genes of 
IncRNAs. Our observations demonstrated that IncRNAs expres- 
sion profiles are likely to provide important insights into 
pathogenesis of HCC. 



Materials and Methods 

Patient sample preparation 

Nineteen parrs of HCC samples and adjacent non-tumor tissues 
were obtained from surgical specimens at General Hospital of the 
People's Liberation Army (Beijing, China) after informed consent. 
Those samples used in this study have been described in previous 
publication [23]. Clear hepatocellular carcinoma was diagnosed 
histopathologicaUy. All these specimens including clear hepatocel- 
lular carcinoma and adjacent non-tumor tissues were snap-frozen 
at liquid nitrogen after excision. Of these samples, 3 pairs (2 
patients infected with Hepatitis B, and one infected with Hepatitis 
C) were used for IncRNA microarray analysis and the remaining 
samples were used for further validation. The clinical character- 
istics of the three HCC patients were shown in Table 1. 

RNA extraction 

Total RNA was extracted from the frozen tissue block using 
TRIzol reagent (Sigma, USA) according to the manufacturer's 
protocol. RNA quantification and quality were assured by Nano 
Drop ND-1000 spectrophotometer. RNA integrity and gDNA 
contamination were tested by denatured agarose gel electropho- 
resis (Figure SI). 

Microarray analysis and computational analysis 

Sample preparation and microarray hybridization were per- 
formed by Kangchen Bio-tech, Shanghai P.R. China. Briefly, 
RNA was purified from 1 |j,g total RNA after removal of rRNA 
(mRNA-ONLY Eukaryotic mRNA Isolation Kit, Epicentre). 
Then, each sample was amplified and transcribed into fluorescent 
cRNA along the entire length of the transcripts without 3' bias 
utilizing a random priming method. The labeled cRNAs were 



Sample Processing 



Expression Profiling 




Bioinfomiatics Analysis 



• Sample preparation 

• Total RNA extraction & QC 



Data collection and Normalization 

Compare for common genes 

LncRNAs classification and subgroup analysis 



GO enrichment analysis 
LncRNA-mRNA co-expression network 



Specific IncRNAs expression 
LncRNAs and their nearby genes 



qRT-PCR Verification 



Figure 1. Flowchart of the experiment. 
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Table 1. Clinical Characteristics of the Patients. 







Patient No. 


Age 


Gender 


Turner Size (cm*cm*cm) 


HBsAg 


HCV-Ab 


Cirrliosis 


HI 


37 


iVl 


1 6*4*8 


Positive 


Negative 


Yes 


H2 


64 


M 


8*5*4 


Negative 


Positive 


No 


H3 


52 


M 


7*7*7 


Positive 


Negative 


Yes 



HBsAg indicates hepatitis B surface antigen; HCV-Ab, hepatitis C virus antibody; IVl, male. 
doi:1 0.1 371 /journal.pone.Ol 01 707.t001 



hybridized onto the Human LncRNA Array v2.0 (8x60 K, 
Arraystar). After having washed the slides, the arrays were scanned 
by the Agilent Scanner G2505B. Agilent Feature Extraction 
software (version 10.7.3.1) was utilized to analyze acquired array 
images. Quantile normalization and subsequent data processing 
were carried out using the GeneSpring GX v 11.5.1 software 
package (Agilent Technologies). DifTerentiaUy expressed LncRNAs 
and mRNAs were identified through fold change filtering (Fold 
Change S2.0 or ^0.5), paired t-test (p<0.05) and multiple 
hypothesis testing (FDR<0.05). P values and FDR were calculated 
by Microsoft Excel and MATLAB respectively. The microarray 
data has been deposited in NCBI Gene Expression Omnibus 
(GEO) and the GEO accession number is GSE55191 (http:// 
www.ncbi.nlm.nih.gov/ geo/query/acc.cgi?acc = GSE55 191). Path- 
way analysis and GO analysis were applied to determine the roles 
of these differentially expressed mRNAs played in these biological 
pathways or GO terms. Differentially regulated mRNAs were 
uploaded into the Database for Annotation, Visualization and 
Integrated Discovery (DAVID, http://david.abcc.ncifcrf gov/) to 
analyze the enrichment of these coding genes. After uploaded 338 
seqnames of differentially regulated mRNAs into Functional 
Annotation Chart, 325 mRNAs were matched to the list. The 
annotation summary results were shown up by this webserver. 

Quantitative real-time polymerase chain reaction 
(qRT-PCR) validation and Statistical analysis 

Total RNA was reverse-transcribed using ImProm II Reverse 
Transcriptase (Promega) according to the manufacturer's protocol. 
Real-time PCR was done with SYBR Premix Ex Taq (TaKaRa) 
on MX3000p instrument according to the manufacturer's 
protocol. Specific primers of each genes are listed in Table SI. 
The relative fold change was calculated using the 2 ^^'^^ method 
normalized to GAPDH. Differences of IncRNAs between tumor 
and normal were analyzed using paired t-test. P value <0.05 was 
considered as statistically significant. 

Co-expression network construction 

The IncRNA-mRNA co-expression network was constructed 
based on the correlation between the differentially expressed 
IncRNAs and mRNAs. The algorithm was quoted from previously 
described report [24]. Red ellipse represents as down-regulated 
IncRNA, green diamond represents as down-regulated mRNA, 
blue ellipse represents as up-regulated IncRNA, black diamond 
represents as up-regulated mRNA. Sohd lines represent positive 
correlation, dotted lines represent negative correlation. 

Results 

Expression profile of IncRNAs and mRNAs in HCC 

Fold-change (Tumor vs Normal) and P value were calculated 
from the normalized expression. Using microarray analysis, we 
identified 214 IncRNAs and 338 mRNAs to be significantly 



differentially expressed in three HCC samples compared with their 
adjacent NT samples (Fold Change &2.0 or £0.5, p<0.05 and 
FDR<0.05). Among these, 166 IncRNAs and 219 mRNAs were 
identified to be consistently upregialated in all three HCC groups, 
and 48 IncRNAs and 1 1 9 mRNAs were consistently downregu- 
lated (Table S2). The expression of ENST00000448647 (Log2 
Fold Change T/N = 7.8561822) was significantly upregulated, 
while the expression of ENST00000504368 (Log2 Fold Change 
T/N= —7.4821358) was dramatically downregulated. The num- 
ber of deregulated IncRNAs and mRNAs varied in different 
patients (Table 2, 3). These IncRNAs may play important roles in 
the occurrence and development of HCC. 

LncRNA Classification and Subgroup Analysis 

According to previous reports, IncRNAs can be classified into 
different subgroups, such as HOX IncRNAs, IncRNAs with 
enhancer-like function (IncRNA-a) and large intergenic iioncoding 
RNAs (lincRNAs) [25,26]. We carried out IncRNA classification 
and subgroup analysis to further investigate potential functions of 
these HCC associated IncRNAs. 

Rinn et al. have identified a total of 407 discrete IncRNAs within 
the four human homeobox transcription factors (HOX) clusters, 
partitioned into 101 known HOX gene exons, 75 introns and 231 
intergenic transcripts [26] . Among them, 6 of coding transcripts 
and 4 of noncoding transcripts were differentially expressed in 
three sets of HCC samples (Table S3). 

Using the GENCODE annotation [27] of the human genome, 
Orom et al. have defined a set of IncRNAs with enhancer-like 
function in human cell lines [25]. Depletion of these IncRNAs led 
to decreasing expression of their neighboring protein-coding 
genes. The profiling data showed that 506 such IncRNAs were 
detected, while 26 of them were differentially expressed. Among 
these 26 IncRNAs, 15 were located in nearby aberrantiy expressed 
coding genes (distance <300 kb) (Table S4, S5), implying that 
these IncRNAs might be regulators of their neighboring coding 
genes in HCC cells. 

By analyzing chromatin-state maps, Rinn and his colleagues 
have identified ~3,300 lincRNAs with clear evolutionary conser- 
vation and association with diverse biological processes, including 
cell proliferation, immune surveillance, ceU-cycle regulation and 
embryonic stem cell pluripotency [28,29]. The profiling data 
indicated that 46 of such lincRNAs were dysregulated in HCC 
samples, and 35 among them were located in nearby aberrantiy 
expressed coding genes (distance <300 kb) (Table S6, S7). 

Bioinformatic analysis of mRNAs expression in HCC 

Up to 12,751 coding transcripts could be detected in three pairs 
of samples through 30,215 coding transcripts probes. 338 mRNAs 
were differentially expressed in all three pairs of samples. Using 
DAVID Functional Annotation Chart, we analyzed the enrich- 
ment of these 338 differentially regulated mRNAs (Table S8). The 
result showed that the most significant functional groups consisted 
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Table 2. Number of aberrantly expressed IncRNAs in microarray 


for three pairs of HCC and 


adjacent non-tumor tissues. 




Sample ID 




Long noncoding RNA 










Fold change 2-4 


Fold change 4-6 


Fold change >6 


Total 


Aberrant IncRNAs 


HI 


up 


1917 


461 


335 


2713 


4369 




down 


1120 


235 


301 


1656 




H2 


up 


1393 


301 


240 


1934 


3312 




down 


1033 


177 


168 


1378 




H3 


up 


1314 


665 


1901 


3880 


7093 




down 


1898 


628 


687 


3213 




doi:l 0.1 371/journal.pone.01 01 707.t002 



of metabolic process, cellular component organization, cell cycle 
and response to chemical stimulus (Figure 2). Liver is a highly 
metabolic organ balancing various biological processes. Its 
disorder is an important cause for the development and 
progression of HCC [30]. Specific nutrients or metabolites are 
associated with increased HCC risk. Recent data have suggested 
that several metabolic regulatory drugs may have the potential to 
decrease the incidence of HCC [31,32]. 

Genomic location of these differentially expressed 
IncRNAs 

The functional characterization of these IncRNAs presents a 
formidable task, and some evidences suggested that IncRNAs 
regualte high order chromosomal dynamics, sub-cellular structural 
organization and telomere biology [33]. LncRNAs have been 
dedicated to originate from complex loci that include IncRNAs 
and associated protein-coding genes. For example, a number of 
IncRNAs have been reported to regulate the expression of 
adjacent protein-coding genes [34,35]. The relationship of 
IncRNAs and its nearby coding gene identified here includes 
bidirectional, exon sense-overlapping, iiitergenic, intro sense- 
overlapping, intronic antisense and natural antiseiise. Among the 
214 consistently regulated IncRNAs in three sets of HCC samples, 
including 29 directional sequences (23 upregulated and 6 
downregulated), 18 exon sense-overlapping sequences (16 upre- 
gulated and 2 downregulated), 120 intergenic sequences (93 
upregulated and 27 downregulated), 4 intro sense-overlapping 
sequences (3 upregulated and 1 down regulated), 29 inronic 
antisense (19 upregulated and 10 downregulated) and 27 natural 
antisense (23 upregulated and 4 downregulated). Such information 
about these IncRNAs and their nearby coding genes might be 
useful to predict their functional roles in HCC (Table S2). 



Construction of the IncRNA-mRNA Co-expression 
Network 

In order to ascertain the correlation between differentially 
expressed liicRNA and mRNA in HCC, we constructed a 
IncRNA-mRNA co-expression network as described [24] . Pearson 
correlation coefficients between all aberrant IncRNAs and mRNAs 
were calculated. Using P-value <0.001 and absolute value of 
correlation coefficient &0.99, we identified 249 pairs of co- 
expressed IncRNAs and mRNAs composed of 131 mRNAs (39% 
of all differentially expressed mRNAs) and 103 IncRNAs (48% of 
all difFerentiaUy expressed IncRNAs), and 146 pairs presented as 
positive correlation (Figure S2). LncRNA-mRNA pairs with the 
highest positive correlation coefficieiicy included BCO 16002- 
ZNF777, ENST00000424412-TUBB2A and CR625009-RPL17. 
ENST00000444546-CES1 and ENST00000433673, uc.472-PLG 
had the highest negative correlation coefficiency. This network 
indicated that one IncRNA could correlate with one to tens of 
mRNAs and so was the mRNA. Increasing evidences indicated 
that IncRNAs play roles in gene expression [36]. For example, 
long intergenic non-coding RNA-p21 activated by p53 mediates 
global gene repression in the p53 response [37]. Thus, we 
proposed that the expression profile of mRNAs and IncRNAs are 
significandy correlated. 

A set of IncRNAs were dysregulated in HCC samples 

To verify the microarray data, we randomly selected fifty 
differentially expressed IncRNAs including twenty up-regulated 
IncRNAs and thirty down-regulated ones and validated their 
expression levels by quantitative RT-PCR (qRT-PCR) in three 
sets of HCC tissues. The result showed that the expression patterns 
of forty IncRNAs were consistent with the microarray data 
(Table 4). Subsequently, eight IncRNAs from these forty were 



Table 3. Number of aberrantly expressed mRNAs in microarray for three pairs of HCC and 


adjacent non-tumor tissues. 




Sample ID 




Message RNA 










Fold change 2-4 


Fold change 4-6 


Fold change >6 


Total 


Changed mRNAs 


HI 


up 


2295 


487 


428 


3210 


6077 




down 


1852 


419 


596 


2867 




H2 


up 


2001 


380 


326 


2707 


5132 




down 


1595 


381 


449 


2425 




H3 


up 


1831 


918 


2402 


5151 


10110 




down 


2709 


885 


1365 


4959 




doi:l 0.1 371 /journal.pone.Ol 01 707.t003 
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A Biological Process 
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Cellular catabolic process 
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C Molecular Function 



Binding 
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Transition metal ion binding 
Oxidoreductase activity 
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Tetrapyrrole binding 
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Counts 



200 
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Figure 2. Gene ontology (GO) enrichment analysis for differentially regulated mRNAs. (A) GO analysis of mRNAs according to biological 
process. (B) GO analysis of mRNAs according to cellular component. (C) GO analysis of mRNAs according to molecular function. 
doi:1 0.1 371/journal.pone.OlOl 707.g002 



evaluated in the nineteen pairs of HCC and adjacent NT samples 
using qRT-PCR (Table S9). The results showed that BC017743, 
ENST00000395084, NR_026591, NR_0 15378 and NR_024284 
were up-regulated and NR_027151, AK056988 and uc003yqb.l 



were down-regulated in HCC samples compared with adjacent 
NT samples (Figure 3). Additionally, microarray data also 
provided the nearby coding gene of these differentially expressed 
bicRNAs (Table SIO). We found that the changing trend of 
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IncRNAs and their nearby coding genes is basically identical 
through qRT-PCR analysis of a selected number of differentially 
expressed IncRNAs and their nearby genes (Figure 4). These 
results provided strong evidence for the point that IncRNAs have 
intrinsic cis-regulatory capacity to their own locus. 

Discussion 

HCC is the fifth most common malignancy in human beings, 
accounting for approximately 90% of primary liver cancers [38]. 
Hepatocarcinogenesis is a complicated biological process charac- 
terized by a myriad spectrum of molecular abnormalities. Over the 
past decades, the molecular mechanism of HCC has been 
extensively investigated. However, the exact pathogenesis of this 
disease is stiU vague. Increasing evidence indicates that IncRNAs 
may play a significant role in regulating gene expression [39]. 
LncRNA expression is de-regulated in many types of cancers, such 
as HOTAIR in breast tumours and metastases, BACEl-AS in 
Alzheimer's disease and Gas5 in breast cancer [15,40,41]. 
Dysregulation of IncRNAs, including H 1 9, HEIH, MVIH, HULC 
and MEGS, has been identified in HCC [16-18,22]. Some of 
these IncRNAs were also identified in our microarry data. For 
example, H19 was downregulated in HI and H2 tissues. 
Therefore, it is a critical step to figure out the expression profile 
of IncRNAs and related mRNAs in HCC in understanding its 
pathogenesis of HCC. 

In this study, we investigated gene expression profiles of HCC 
by using IncRNA microarray. Compared with previous report 
[18], which showed hierarchical clustering analysis of 254 mRNAs 
and 174 IncRNAs that were differentially expressed between five 
pairs of HCC samples and nontumor samples (all samples using 
here are HBV-related), we analyzed three pairs of HCC and 
adjacent non-tumor samples (two patients infected with Hepatitis 
B, and one infected with Hepatitis C) and identified 214 IncRNAs 
and 338 mRNAs abnormally expressed in all three HCC tissues 
(Fold Change >2.0, P<0.05 and FDR<0.05) which is partly 
different from previous results. This is due to differences in 
cancerous tissues. Moreover, the microarray which we used is 
human IncRNA microarray version 2.0, whereas the microarray 
which was used in the above report is human IncRNA microarray 
version 1.0. Compared with version 1.0, version 2.0 contains more 
comprehensive and reliable array content, most extensive and 
updated coverage available, specific exon or splice junction 
probes, efficient and robust labeling system and systematic 
IncRNA classification. According to the absolute expression profile 
in three pairs of sample, the numbers of significantly differentially 
expressed IncRNAs and mRNAs (Fold Change £2.0 or <0.5) are 
624 and 1050 separately. After statistical analysis (the threshold for 
deregulated genes was p<0.05 and FDR<0.05), this number 
decreased to only 214 and 338. This phenomenon might be 
because the gene expression level varies considerably in different 
patients. Therefore, this statistical analysis will lead to the results 
that many important differentially expressed genes have been 
excluded. 

Yang F et al. have confirmed that AY129027, uc002pyc, and 
DQ786243 were overexpressed, whereas the expression of 
AK055007 and AK123790 was decreased in 50 paired HCC 
compared with adjacent NT samples using RT-PCR [18]. In our 
study, we demonstrated that BC017743, ENST00000395084, 
NR_026591, NR_0 15378 and NR_024284 were up-regulated, 
whereas NR_027151, AK056988 and uc003yqb.l were down- 
regulated in nineteen pairs of HCC samples compared with 
adjacent NT samples using qRT-PCR. Moreover, expression of 
seven IncRNAs was significantly correlated with their nearby 
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Figure 3. Distribution of eight incRNAs in nineteen pairs of HCC samples. Selected eight differentially expressed IncRNAs were validated by 
qPCR in 1 9 pairs of HCC samples. The heights of the columns in the chart represent the relative concentration of each IncRNA validated across the 1 9 
patients. 

doi:l 0.1 371/journal.pone.01 01 707.g003 



coding genes which have not been described before. For example, 
NR_026591 is a 452^ bp exon sense-overlapping IncRNA 
transcribed from chr5. It is one of the alternatively spliced 
transcript variants of the RAD 1 gene, which encodes a component 
of the 9-1-1 complex, a heterotrimeric cell cycle checkpoint 
complex involved in stopping cell cycle progression. NR_026591 
and BRIXl make up a bidirectional transcript pair, expression of 
which are initiated in opposite directions and in close proximity 
[42]. QRT-PCR results showed that the expression trend of 
protein-coding genes (RADl, BRIXl) and IncRNA (NR_026591) 
was consistent. We speculate that NR_026591 may have cis- 
regulatory capacity to their own locus. All the seven IncRNAs 
(Table SIO) are uncharacteristic and associated biological func- 
tions of them are unclear. The present data will be beneficial to 
more in-depth exploration of HCC. 

LncRNAs are involved in transcriptional regulation of many 
protein-coding genes. In order to obtain novel insight into the 
function of IncRNAs, we further combined differentially expressed 
IncRNAs with differentially expressed mRNAs to construct a co- 
expression network. Owing to differences of abnormally expressed 
IncRNAs and mRNAs, the kicRNA-mRNA co-expresion network 



is different from the network described in the previous report. We 
found that many IncRNAs were significantly correlated with 
multiple protein-coding genes. For example, this network indicat- 
ed that ZNF777 was associated with BC016002. BC()16()()2 is a 
960 bp intergenic IncRNA located on chr7. ZNF777 is involved in 
transcriptional regulation due to its DNA binding function which 
is the most significant functional group in the result of GO analysis 
of differentially regulated mRNAs. Therefore, we hypothesized 
that ZNF777 may be the direct or indirect target genes of 
BCO 16002. This network also indicated that CESl are predicted 
to be associated with ENST00000444546. ENST00000444546 is 
a 1150 bp intergenic IncRNA located on chrl. CESl encodes a 
member of the carboxylesterase large family which are responsible 
for the hydrolysis or transesterification of various xenobiotics. This 
enzyme is the major liver enzyme and functions in liver drug 
clearance. Drug clearance which is an important characteristic of 
CESl may be critical for the function of ENST00000444546. 
ENST00000433673 is a 2690 bp intergenic IncRNAs located on 
chr2 and uc.472 is a 202 bp exon sense-overlapping IncRNA 
located on chrX. Their correlated gene in co-expression network is 
PLG. Protein encodes by this gene is plasminogen which is 
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Figure 4. Consistency of IncRNAs and tKieir nearby genes expression level. Selected seven difFerentially expressed IncRNAs and their nearby 
genes were validated by qPCR in 1 9 pairs of HCC samples. The heights of the columns in the chart represent the log2 fold changes (T/N) in expression 
across the 19 patients for each of the seven IncRNAs validated. 
doi:l 0.1 371/journal.pone.01 01 707.g004 



activated by proteolysis and converts to plasmin and angiostatin. 
Plasmin is an essential protease in a variety of cellular processes 
including embryonic development, tissue remodeling, tumor 
invasion and inflammation. Angiostatin is an angiogenesis 



inhibitor that blocks neovascularization and growth of primary 
and metastatic tumors. This characteristic indicates that 
ENST00000433673 and uc.472 may be associate with PLG to 
exert its function in HCC. The IncRNA and mRNA co-expression 
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network provides a strong foundation for predicting the function of 
IncRNA. The malfunction of regulation in this co-expression 
network may be a possible step for development and progression of 
HCC. The precise molecular mechanisms have yet to be defined. 

In conclusion, we described the expression profile of IncRNAs 
and mRNAs in HCC by microarray and identified 214 IncRNAs 
and 338 mRNAs abnormally expressed in all three HCC tissues 
(Fold Change >2.0, P<0.05 and FDR<0.05). After qRT-PCR 
verification, we demonstrated for the first time that BC017743, 
ENST00000395084, NR_026591, NR_015378 and NR_024284 
were up-regulated, whereas NR_027151, AK056988 and 
uc003yqb.l were down-regulated in nineteen pairs of HCC 
samples compared with adjacent NT samples. Moreover, expres- 
sion of seven IncRNAs was significantiy correlated to their nearby 
coding genes, suggesting these IncRNAs may play roles in HCC 
cells by regulating their nearby coding genes. For these 
dysregulated IncRNAs, we carried out IncRNA classification and 
subgroup analysis, and constructed IncRNA-mRNA co-expression 
network. Such information would be used for investigating the 
functions of these IncRNA in the occurrence and development of 
HCC. Our results revealed a penetrating layer of IncRNA of 
biological significance in pathologies of HCC, thereby pointing 
intriguing directions for further research. 
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